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Q . ABSTRACT 

o : 

^ . Gamma-ray bursts are promising tools for tracing the formation of high red- 

! shift stars, including the first generation. At very high redshifts the reverse shock 

^ ! emission lasts longer in the observer frame, and its importance for detection and 

OO '. analysis purposes relative to the forward shock increases. We consider two differ- 

^ ! ent models for the GRB environment, based on current ideas about the redshift 

O ! dependence of gas properties in galaxies and primordial star formation. We cal- 

Q ! culate the observed flux as a function of the redshift and observer time for typical 

r-| ! GRB afterglows, taking into account intergalactic photoionization and Lyman-a 

^! absorption opacity as well as extinction by the Milky Way Galaxy. The fluxes 

O I in the X-ray and near IR bands are compared with the sensitivity of different 

to '. detectors such as Chandra, XMM, Swift XRT and JWST. Using standard as- 

i sumptions, we find that Chandra, XMM and Swift XRT can potentially detect 

.!-H i GRBs in the X-ray band out to very high redshifts z > 30. In the K and M 

^ i bands, the JWST and ground-based telescopes are potentially able to detect 

I GRBs even one day after the trigger out to 2; ~ 16 and 33, if present. While the 

X-ray band is insensitive to the external density and to reverse shocks, the near 
IR bands provide a sensitive tool for diagnosing both the environment and the 
reverse shock component. 

Subject headings: GRB: massive stars; Cosmology: high-redshift; Spectroscopy: 
x-ray-NIR 



1. Introduction 



Gamma-ray bursts are thought to be associated with the formation of massive stars 
(van Paradijs, Kouveliotou, & Wijers 2000). The evidence for this has been mainly in the 
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class of long bursts, of 7-ray durations in excess of 2 seconds, making up two-thirds of the 
GRB population, which are the only ones so far for which X-ray, optical, IR and radio 
afterglows, as well as redshifts, have been measured. The strongest evidence yet comes from 
the recently confirmed association of long GRBs with corc-collapsc supernovae (Stanek et 
al. 2003; Hjorth et al. 2003; Uemura et al. 2003; Price ct al. 2003). Short bursts, of 
durations less than 2 seconds, even if produced e.g. by neutron star mergers, would similarly 
be associated with massive star formation, and one expects the rate of occurrence of GRBs 
with redshift to follow closely the massive star formation rate. In currently favored AGDM 
cosmologies, star formation should start at redshifts higher than those where proto-galaxies 
and massive black holes at their centers develop (Miralda-Escude 2003). Thus, GRBs could 
trace the pre-galactic star formation era preceding quasars. 

Recent cosmic microwave background anisotropy data collected by WMAP reveal that 
the first objects in the Universe should be formed around z ^ 18 (Bennett et al. 2003). This 
is consistent with the theoretical modeling of the first star formation (Abel et al. 1998, 2000, 
2002; Bromm et al. 1999). There is also indirect observational evidence for high-z GRBs. 
E.g., empirical relations have been found between the GRB luminosities and other measured 
quantities, such as the variability of the gamma-ray light curves (Fenimore & Ramirez-Ruiz 
2000) and spectral lags (Norris et al. 2000). By extrapolating these empirical laws to a 
larger burst sample (e.g. the BATSE data), it is found that many BATSE bursts would be 
expected to have 2; > 6 (Fenimore & Ramirez- Ruiz 2000). 

The discovery of the highest redshift quasars, such as the current record holder at 
z — 6.43 (Fan ct al. 2003), grows increasingly difficult because the quasar formation rate 
drops rapidly at higher redshifts, peaking between redshift 2 and 3. Very few galaxies can 
be seen above 2; > 6, which is also consistent with the upper limit for the redshift of galaxy 
formation Zgai < 9 based on theoretical analysis (e.g. Padmanadhan 2001). Although young 
galaxies may exist at very high redshifts, they are likely to be too faint to obtain good spectra 
(Haiman & Loeb 1997). On the other hand, the extreme brightness of GRBs during their first 
day or so make them the most luminous astrophysical objects in the Universe. Thus, GRBs 
appear to be promising tools to explore the very high redshift Universe (Miralda-Escude, 
1998). 

The natural question which needs to be quantified is the degree of detectability of GRBs 
with current or future detectors, if they occur at much higher redshifts than those currently 
sampled. Lamb & Reichart (2000) used specific templates such as GRB 970228 observed 
at one day to estimate the highest redshifts at which such bursts could be observed using 
Swift. Ciardi & Loeb (2000) calculated the flux evolution with redshift of common GRBs and 
discussed the flux change with redshift at several epochs in the infrared bands. These papers 
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considered only forward shock radiation as known before 2000 and some effects of the galactic 
mean density evolution but did not consider the primeval star-formation environment. 

In this paper we have calculated the flux evolution of typical GRBs based on current 
knowledge about GRB physics in a more realistic way. Among the refinements introduced 
are: (1) The contribution from reverse shocks is considered as a crucial element. This 
should be very important for the early afterglow in the rest-frame (which at high redshifts 
gets dilated to longer observed times). Therefore, we expect that at higher redshifts the 
possibility of observing the reverse shock is much increased. (2) We have taken up to date 
GRB parameters, e.g. incorporating new estimates of the typical magnetic equipartition 
parameter about one order or more magnitude smaller than the electron parameter tf. in 
the forward shock, and a possibly higher in the reverse shock. This has a significant effect 
on the GRB evolution. (3) We consider GRB external densities motivated both by views 
on the typical protogalaxy density evolution with redshift, and by views on the conditions 
around the first stars to form in the universe in the pre-galactic era. (4) We consider both 
the Lyman-a and photoionization absorption as well as our own galactic extinction. (5) We 
compare the expected fluxes in the X-ray and near IR bands to the sensitivity of various 
detectors such as Chandra, XMM, Swift XRT and JWST. 

In §§2.1 and 2.2 we outline the basic forward and reverse shock flux calculations, the 
details of which are given in an appendix. We discuss the GRB density environment in §2.3, 
and the intergalactic and galactic absorption effects are estimated in §2.4. In §§3.1 and 3.2 
we discuss the optical/IR and X-ray flux dependence on redshift, respectively, at various 
observer times, including the dependence on external density. We compare these to the 
Swift XRT, Ghandra, XMM and JWST sensitivities for the detection of GRBs at different 
redshifts. We summarize the numerical results and discuss the implications in §4. 

2. Afterglow chciracteristics 

2.1. Forward Shock 

We assume that the shock-accelerated electrons have a power-law distribution of Lorentz 
factors 7e with a minimum Lorentz factor 7^: N{'^e)dle oc ^^d^y^ile > 7m- We also define 
a critical Lorentz factor 7c above which the electrons cool radiatively on a time shorter 
than the expansion time scale (Meszaros, Rees, & Wijers 1998). This leads to the standard 
(forward shock) broken power law spectrum of GRBs (Sari, Piran, & Narayan 1998). In the 
fast-cooling regime, when 7^ > 7c, ah the electrons cool rapidly down to a Lorentz factor 
~ 7c and the observed flux at frequency v is 
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Fu = i^.,m,/ < {y/ycjY^'^ ucj <y< l^mj (1) 

Hereafter the subscripts 'f and 'r' indicate forward and reverse shock, respectively. 

In the slow-coohng regime, when 7c > ^rm only electrons with 7e > 7c cool efficiently, 
and the observed flux is 

where Fy^^nj is the observed peak flux at the observed frequency v = min(z/cj, i^mj), while 
u^j and i/c,f are the observed frequencies corresponding to 7^ and 7c, respectively. Syn- 
chrotron self- absorption can also cause an additional break at very low frequencies, typically 
about < 5 GHz, in the radio range. Since here we focus on the IR and X-ray ranges, we will 
not consider this low- frequency regime in our calculations. For a fully adiabatic shock, the 
evolution of the typical frequency and peak flux are given by Sari et al. (1998): 

i/cj = 2.5 X 10^2^-'f Eg-'/'n-H"'/' X (1 + 

urnj = 4.6xl0^^e;/>^Ejft-^/^x(l + z)V^Hz, (4) 
Fy,m,f = l.lxWh]^^fE52n'/^D2s{z)-'' x{l + z) fxJy . (5) 

Hereafter the quantities without the subscript 's' are in the observer frame, and the quantities 
with the subscript 's' are for the observer in the local frame of the source, which is connected 
with the observer frame quantities with a certain power of (1 -|- 2;). The source is assumed 
at a luminosity distance Dl{z) — 10^^D28(-s) cm, and es and eg are the fraction of the shock 
energy converted into energy of magnetic flelds and accelerated electrons, respectively. The 
time is taken in units of t = lO^ts s (~ 1 day), £'52 = ergs is the isotropic equivalent 

energy of the GRBs, and n is the particle density in units of cm~^ in the ambient medium 
around the GRB. 



2.2. Reverse Shock 

As GRBs are measured at increasingly larger redshifts, a given constant observer time 
corresponds to increasingly shorter source frame times. This is favorable for observing at 
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very high redshifts the evolution of phenomena which happen only in the earliest stages 
of the GRB, such as the reverse shock emission. So far, the reverse shock emission has 
been observed in only three GRBs in the optical band: GRB 990123 (Akerlof et al. 1999), 
GRB 021004 (Fox et al. 2003a) and GRB 021211 (Fox et al. 2003b; Li et al. 2003). 
At these early epochs, the reverse shock emission makes a significant contribution to the 
overall flux of the GRB afterglow. A description of the reverse shock spectrum is however 
more complicated than that of the forward shock. It depends on two factors: (1) whether 
one is in the thick shell or thin shell case, and (2) the ratio of the crossing time of the 
reverse shock across the shell to the observing time. We consider a rclativistic shell with 
an isotropic equivalent energy E and initial Lorenz factor r] = L^/Mc^ expanding into 
a homogeneous interstellar medium of particle number density n. In the local frame, we 
can define a deceleration timescale when the accumulated ISM mass is l/r] of the ejecta 
mass, tdec,s = \{'iE / AurfnufipC^Y^'^ /2rfc\^ which is the conventional deceleration timescale. 
A critical initial Lorenz factor rjc can be defined by the condition that the deceleration time 
tdec,s equals the intrinsic (i.e. central engine dominated) duration T, of the gamma-ray burst, 
which is - '22%.&Eli^n-^I^T;l'^ . The thick shell case occurs when the duration Tg > tdec s-, 
and the thin shell case occurs when Tg < tdec,s- The time taken by the reverse shock to cross 
the shell is defined as tx,s — T^^^{tdec,s,Ts)- In the observer frame, tx — tx,s(l + z). For 
observation times t < tx- the reverse shock emission spectrum qualitatively resembles the 
forward shock spectrum. However, for t > tx, we take the approximation that there is no 
reverse shock emission above z^cr, since all electrons have cooled below that energy. Thus 
the reverse shock spectrum at t > tx is, in the fast cooling case 

{{l^/l^Cry^^ ^ < ^c,r 

Ucj- <1^< l^m,r (6) 

l^m,r < 

In the slow cooling case the reverse shock spectrum is 

ij^/j^m,r)~^'^^^^ J^m,r < < ^cr (7) 

Vc,r < y 

where i^cr, i^m,r and Fi,^jn,r refer here to the reverse shock cooling frequency, typical frequency 
and peak fiux, respectively. Since these quantities are usually different in the reverse and in 
the forward shocks, and have different functional forms and time evolution dependences in 
the thick and thin shell cases, the specific shock and shell cases will be differentiated in the 
treatment below. 

Kobayashi (2000) has given expressions of cooling frequency, typical frequency of elec- 
trons and peak fiux in the reverse shock. Motivated by recent observations of prompt fiashes. 
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a set of relations linking the cooling frequency, typical frequency and peak flux in the reverse 
and forward shocks at the crossing time was proposed by Kobayashi & Zhang (2003) and 
Zhang, Kobayashi, & Meszaros (2003). The flux calculated with these two different sets of 
formulae are consistent within a 10% error range. Here we use the relationships as discussed 
in the last two quoted references, 

l^m,r{tx) / 2 / \-2-T-> ^c,r\tx) ^ _3 -^!^,m,r(^x) / 2 / \t> fo\ 

- — TTT = iixh) 7^B , , . = 7^B , — 7-^ = [lxh)T^B (8) 

where 

7x = min(r7, rfc) , and TIb = Br/Bf = {eB^r/eBjf^ ■ (9) 

Here TZb reflects a possible stronger B fleld in reverse shock, as inferred from the analyses 
of the GRB 990123 and GRB 021200 data (Zhang et al. 2003). In our calculations we set 
TZb = 1 as the standard case, and take TZq = 5 as an alternative option. As an example, 
when the observer time is larger than the crossing time, t > ix=iiiax(T, i^ec), i-e. the fast 
cooling case, the observed cooling frequency, typical frequency and peak flux of the reverse 
shock are 

^c,r = {txltf'^nB-^V,j{t^) = 2.5 X 1012,^3/2^-1/2^-1^-3/2^^^^ ^ ^ zY^'^Ub-^ Hz, 
l^m,r = 7x V(^x A)'/'7^BZ/„J(^x) = 4.6 X lO^^" V^^Jee^sf ^5 X (1 + ^Y^'^^B Hz, (10) 

F,,m,r = ^lr\t^/t)nBF,^„,j{t^) = 1.1 X lQh'ilE^2n^'^D2s{z)-^ X (1 + z)^lr]-H^,^tl^nB i^Jy. 

In the appendix we give further details of the expressions for the flux evolution of the 
forward and reverse shocks in the thin and thick shell as well as in the fast or slow cooling 
cases. 



2.3. GRB density environment 

The typical environments considered for GRBs arc either the (approximately) constant 
number density case uq ~ constant (i.e. independent of the distance r from the center for 
the burst), or a power law dependence as might be expected in the stellar wind from the 
progenitor, e.g. n oc r"^ (Meszaros ct al. 1998; Dai & Lu 1998; Chevalier & Li 1999; 
Whalen, Abel, & Norman, 2003). In our calculation, for simplicity we consider only the 
first case of n ~ constant, which appears to satisfy most of the observed cases that have 
been analyzed (Panaitescu & Kumar 2001,2002; Frail et al. 2001). While this density is 
different for different bursts, we can assume a typical average value no for n at redshift 
z — 1. One has to consider then how this typical density might evolve with redshift. We 
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concentrate on two very different types of dependencies, motivated by different physics. (1) 
Based on hierarchical models of galaxy formation (Kauffmann, White, & Guiderdoni 1993; 
Mo, Mao, & White 1998), the mass and size of galactic disks is expected to evolve with 
redshift (Barkana & Loeb 2000). For a fixed host galaxy mass, this yields n{z) = no(l + z)^ 
(Ciardi & Locb 2000). (2) Recent numerical simulations of primordial star formation indicate 
that the particle number density around the first stars at very high redshift could be in the 
range 1 Uq < 10~^ cm~^ (Whalen et al. 2003), approximately independent of redshift 
because of strong radiation pressure from the central massive star, which dominates and 
smooths any variations in the original galactic number density around the stars. The size 
scale of this region of dominance is about several parsecs, which is > the length scale of 
typical afterglows. Here we assume that, for this case (2), this stellar dominance applies to 
all GRBs originating from massive stars, so the number density in the relevant region around 
the GRB is the same constant at all redshifts, i.e. n{z) = n. Thus, the two density cases 
considered are 

, , _ f (1 + 2;)° constant density model , , 

niz) = Hq Xs. . . (11) 

1^ [1 + z) density evolution model 

Here Uq is normalized by no = 1 cm~^ at 2; = 1, noting that uncertainties in the primordial 
star calculations could make this as low as 10~^ cm~^. This number density n refers to the 
local ISM density in the immediate neighborhood of the burst. 



2.4. Intergalactic and galactic absorption 

As it propagates through the intergalactic medium (IGM), the afterglow radiation from 
a burst occurring at some redshift z is subject to several absorption processes. The most im- 
portant are Lyman-a absorption, photoionization of neutral hydrogen, and photoionization 
of He H. At very high redshifts, before the intergalactic medium becomes re-ionized, which 
may be taken to occur between the limits Zi > 6.3 (Fan et al. 2001; Miralda-Escude 2003; 
Onken & Miralda-Escudc 2003) and Zi ~ 17 ± 5 (Spergel et al. 2003), most of the mass as 
well as most of the volume of the IGM is in the form of neutral gas. At redshifts below this, 
after rc-ionization by the first stars or galaxies, an increasing fraction of the IGM volume 
becomes ionized, interspersed with clouds of neutral gas associated with the halos of proto- 
galaxies, which continue to absorb radiation. The exact distribution of clouds as a function 
of redshift is not well known, but estimates of the effective number are obtained by counting 
the numbers of absorption line systems in quasar spectra. These are used for calculating 
the effective absorption optical depth at redshifts below the reionization redshift. Below the 
rcionization redshift, the photoionization opacity by HI is given by Madau, Haardt, & Rees 
(1999), based on the observed absorber distribution in the spectra of high- redshift quasars. 
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The Lyman-o; absorption optical depth can be obtained in a similar way. Above the reion- 
ization redshift, both the photoionization and Lyman-a opacities are obtained by means of 
an integration through the neutral gas between the reionization redshift and the redshift at 
which the GRB is located (Barkana & Loeb 2001). 

At high redshifts, intergalactic He II becomes important at rest-frame energies > 54.4 
eV, where the effects of hydrogen photoionization are still important (Perna & Loeb 1998). 
However, the combined effect of the cross sections and the abundances, as well as the hardness 
of the ionizing spectra combine together to make He 11 the dominant opacity at observed 
photon energies hu > 54.4eV/(l + z) for sources located at 2; > 3 (Miralda-Escudc 2001). 
Bluewards of this energy, as the cross section drops as u^^, He 11 photoionization is the 
last process to become optically thin, and is therefore the dominant IGM constituent which 
determines the re-emergence of the source spectrum at frequencies above the blue end of the 
Gunn-Peterson trough. Adopting current values of the cosmological parameters, this occurs 
(Meszaros & Rees 2003) at soft X-ray energies of huf ~ 0.2 keV or i/^ ~ 5 x 10^^ Hz. 

Absorption by our own galaxy also becomes important in the UV and soft X-ray band. 
The combined cross section including galactic metals aph is given by Morrison & McCammon 
(1983). The optical depth is given by r = (TphNn, Galaxy, where NH,Gaiaxy is the equivalent 
column density along the line of sight, which varies depending on the galactic latitude. Here 
we set the column density to be 2 x 10^° cm^^, typical of moderately high latitudes, which 
becomes optically thin at energies > 0.2 keV, comparable to the effects discussed above for 
the intergalactic He II. 

Thus, one expects that between the Lyman-a frequency corresponding to the source 
frame and approximately 5 x 10^^ Hz (below which the galactic extinction for the above 
column density becomes large), the flux observed from a high redshift GRB will be totally 
suppressed. Outside this range, the observed flux is much less affected by the intergalactic 
and galactic absorption. 

3. Initial Conditions and Numerical Results 

In our calculations, the nominal GRB parameter values adopted are an isotropic-equivalent 
energy E^2 = 10, shock parameters eg = 0.1, e^^,. = 0.025, the ratio of magnetic flcid strength 
in reverse and forward shocks T^g = B^/Bf = 1 or 5 (e^j = 0.025 or e^j = 0.001, respec- 
tively), and an initial Lorentz factor r] — 120. The GRB duration is assumed to be = 10 s 
in the source frame. The deceleration time tdec,s = [{^E / ATxrfnnipC^Y/'^ /2rfc\ in the source 
frame is determined mainly by GRB intrinsic parameters, except for the external ISM den- 
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sity n, which can depend on redshift in one of the scenarios considered. Substituting the 
parameters for tdec,s, we have tdec,s = 119.6n~^/^seconds. Therefore, for the n=const scenario, 
the reverse shock is exclusively in the thin shell case; for n oc (1 + 2;)^ scenario the reverse 
shock will be in the thin shell case below some redshift, and above that redshift it will be in 
the thick shell case. We take a specific case where the reionization redshift of the Universe 
is at Zi = 15, compatible with the WMAP value of Spergel et al (2003). As examples, we 
considered the burst properties at various observer's times, e.g. 10 minutes, 2 hours and 1 
day. The results are presented in Figs. 2 and 3, discussed below. 

3.1. Light Curve 

To show the distinct effects of the reverse and the forward shock on the flux behavior, 
we show the light curves for two different observational bands, V and K, here taken at a 
nominal redshift z — 1. (Figure 1). As can be seen, the light curve evolution can be divided 
into 3 stages: (1) Forward shock dominant, before the reverse shock peaks. The timescale for 
this stage is relatively short. (2) Reverse shock dominant, after the reverse shock emission, 
which increases very quickly, exceeds the forward shock emission. The flux peaks when the 
reverse crosses the fireball shell. (3) Forward shock dominant. In this stage, after the reverse 
shock peak, the reverse shock emission decays very quickly and falls below the forward shock, 
so the forward shock is again dominant. 

3.2. Infrared Flux Redshift Dependence 

To test the self consistency of the code with the present optical observation like ROSTE, 
we plot the light curves in V band (Figures 2(a) and 3(a)). On these two figures the dashed an 
solid lines correspond to the sensitivities of ROTSE at very early and late times, respectively. 
The light curves show that we don't expect to see many optical flashes, especially for high-^ 
GRBs. This is consistent with the present observations. 

From the flux evolution equations (see appendix), it is seen that in the regime where 
the observing frequency is above the cooling frequency, u > u^j, the observed flux for the 
forward shock component is independent of the ISM number density (see (A4) & (A5)). We 
can define a critical redshift Zc, such that for z > Zc the GRB afterglows for the forward 
shock component are in the density-independent regime (see eq.(3)). 

(l + Zc) = 3Axl0'^{eB,f/0M)-^E^2^n-h.^\\/2.2fj,mf 

= 1.1 X 10^(eB,f/0M)-^E^^^n-H^\leV/hu)-^ (12) 
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Prom the above equation (12), we see that the dependence of Zc on eg,/ is very sensitive, 
(X e^j. If we set e^j = 0.1, we obtain the equations given by Ciardi & Loeb (2000). Taking 
e^j smaller, the redshift Zc can increase substantially. That is one of the main reasons why 
our curve of flux vs. redshift differs from that of Ciardi and Loeb (2000). 

For the description of reverse shocks there are four relevant cases, depending on whether 
one is in the thin or thick shell limit, and on whether the times considered are before or after 
the shock crossing time. However, the cooling frequency evolution can be approximated by 
z/c,r oc t^'^/^ for observer times t > tx = max{tdec, T). Also, if the observed frequency is larger 
than the cooling frequency v > Vc,r-, the reverse shock emission disappears. Hence, we can 
define another critical redshift z^. at which the reverse shock emission disappears, 

(1 + Zr) = 3.33 X lQ-^{eB,r/^-^l)~^E^^n-'^{t^/m)Hf{\/2.2iim)'^nB~^ (13) 

which is a lower limit for n oc (1 + zY and is an upper limit for n = no =constant, i.e. for 
z > Zr when n (x {1 + z)^ or for z < Zr when n=const, there is no reverse shock emission 
(see below). 

The two critical redshifts can be connected by the relations u = Vc,r = (^x/^)^''^^e~^^'c,/(-2r) 
and u = Ucjizc). Cancelling out u and substituting the expression for Ucj, we obtain the 
relation 

{l + Zr) = {tx/tfnB-'{l + Z,), (14) 

and we have the inequality 

Zr < Zc (15) 

since TZb > 1, and t > tx by default, i.e. Zr is defined for t > tx- 

Using the parameters above and an observer frequency u = 1.36 x 10^^ Hz (6.3 x 10^^ 
Hz) or A = 2.2 /xm (4.8 /xm), corresponding to the K-band (M-band) at observer times t — 10 
mins, 2 hour and 1 day, we have {TZjs = 1) 

^ ^ r 2.2, 0.1, (3.0, 0.4, 0) for n cx (1 + z)^ ; 

~ \ 0, 120.6, oo (0, 24.5, oo) for n = 1 cm=^ . ^ ' 

From equation (11) we have, for the n oc (1 + z)^ case, 

l^c,r,evolv = 2.5 X 10'%f E-^'^\l,^^^^i^tf%,c,s,5 X (1 + ^) -^/^^b'^Hz OC (l + ^)-=^9/^ (17) 

and for the n — const case 

l^c,r,const = 2.5 X lO'^-J^^' E-,'/\^,^^^J,'%,,,,,, X (1 + zY/'TZs-'Rz OC (1 + z^/^ (18) 
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Prom equations (13), (17), (18), we find some interesting differences between the two 
density profile cases. We discussed that for i/ > u^r there is no emission from the reverse 
shock. However, the above behavior of Ty'c,r;evoiv (1 + z)~'^^^^ for n oc (1 + z)^ and h'c,r,const oc 
(1 + z)^/"^ for the n = Uq = constant case has some other consequences. For the same burst 
parameters, if there is no reverse shock emission at some observer time t for z = 0, in the 
n oc (1 + z)'^ case this implies that there will be no reverse shock emission at this same 
observer time at any redshift. For the n — constant case, however, the chances are that the 
reverse shock emission will be observable above some redshift, because iy'c,r,const increases with 
redshift. This is caused by the effect of time dilation increasing as the redshift increases, 
which means that the same observer-frame time corresponds to earlier and earlier source- 
frame times. In the n oc (1 -|- z^ case, Ucr.evoiv oc (1 -|- z)"^^/^, and we can expect that z/c,r 
will decay with z quickly below z/ even if z/c,r is much larger than u at low redshift. So, 
in this case, we can only observe the reverse shock emission at relatively low redshifts. On 
the other hand, in the n = constant case we can observe the reverse shock emission at all 
redshifts if there is emission at low redshifts. We can see this from the fluxes in figure 2. If 
we substitute the values for the relevant parameters, we get z < Zr — 2.1 and z < Zr — 
at i = 10 mins and t — 1 hour respectively for the n oc (1 -|- z)^ case, whereas for the 
n — const case, z > Zr — and z > Zr — 1975 for those two corresponding times, and the 
emission from the reverse shock is observable in K band. This property provides one way of 
distinguishing these two different density profile regimes, based on the redshift distribution 
of the occurrence or absence of a reverse shock component. 

From equations (14,15) we see that iit > for the n oc {l + z^ case, the reverse shock 
emission is absent already at redshifts lower than those beyond which the GRB emission 
would be in the density- independent regime. On the contrary, for the n=const case, reverse 
shock emission exists when the GRB emission is in the density-independent regime. Using 
this characteristic, we can again constrain the density profile around the GRBs. 

Looking at Figure 2 which is the standard case of TZq = 1 in our calculation, several 
features can be noted: (1) At early times, e.g. t = 10 mins and t = 2 hour, we can 
differentiate the constant density profile from the evolving density profile in both K and M 
bands. However, at late times it becomes difficult to do so in both bands, although the total 
flux in M band are somewhat different for both density profiles at relatively low redshifts. 
(2) For early observer times, e.g. t — 10 mins, the amplitude of the total flux in the evolving 
density n oc {1 + z)'^ case at low redshifts shows some pronounced and complicated changes 
with redshift, as opposed to a more monotonous behavior in the constant density case. The 
changes in the former are caused by the transitions from one regime to another by the forward 
shock. Because in the evolving density example reverse shock emission disappearance above 
some very low redshift, over the redshift considered here forward shock emission component 
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is dominant. On the other hand, in the constant density case, over the entire redshift range 
considered reverse shock emission is dominant. (3) The break in the Ught curve for the n = 
constant case at t = 10 minutes is caused by the transition from t > tdec,s to t < tdec,s- (4) 
There is a sharp dechne in the emitted flux in hght curves at redshift 2; ~ 17 for K band 
and at z ~ 36, which are caused by the lyman-a and photoionization absorption of neutral 
hydrogen in IGM. 

We also considered the 7^^ = 5 case (Figure 3) which indicates that magnetic field in 
the reverse shock is much stronger than that in the forward shock. Other parameters in this 
case arc same as those in Figure 2. The most distinct feature for this case from the standard 
case is that there is one jump around redshift ^ ~ 2 in the evolving density case. This jump 
is caused by the disappearance of reverse shock above some redshift. Based on the equations 
(A12-A17), the flux ratio between reverse and forward shock increases with an increasing 
TZb- So when the reverse shock emission disappears, we can expect a sudden jump in the 
light curve as a function of redshift. 

We further tested a different normahzation, i.e. uq — 0.01 cm~^, still for TZb — 5. 
The same total forward plus reverse shock flux for the two density profiles for this case are 
shown in Figure 4. An obvious feature is that the reverse shock is much more prominent 

for the lower normalization density (no = 0.01 cm~^) case than for the higher normalization 
density (no = 1 cm~^) case, in the density evolution model. From equations (A4, A5, A20 & 
A19), for slow cooling, the ratio of reverse shock to forward shock emission is proportional 
to n'^^/^^^/'^^ = n^^/^ or n^^/'^^'^/'^) = n^^/^. Therefore, the reverse shock emission becomes 
more prominent for a decreasing number density. However, we also note that the reverse 
shock emission is smaller than the forward shock fiux in some cases. This is because the 
reverse shock is in the i'm,r < ^ < ^c,r regime before the crossing time. So the ratio of the 
reverse shock to forward shock fiux is proportional to n^^^+^^Z^. When the number density 
is smaller than 1 cm~^, we can expect the reverse shock fiux to be smaller. The same is also 
expected for the constant density case at an early observer time. So in this case these two 
density models can be easily differentiated from each other. 



3.3. X-ray Flux Redshift Dependence 

The X-ray band fiux evolution and its redshift dependence is simpler than in the 0/IR 
bands because the reverse shock emission is generally negligible, and we need only consider 
the forward shock emission. One obvious characteristic of Figure 5 is that the fiux from 
the two different density profiles are the same at all the redshifts for a given time, because 
the emission in both cases is in the density-independent regime, being above the cooling 
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frequency Vcj. Based on equation (12), for an arbitrary choice of low X-ray energy of 0.1 
keV and using the other parameters above, we obtain a critical redshift Zc = 0.83 (n=const) 
or Zc = 0.98 (n oc (1 + z)^) where the GRBs change from the density-dependent to the 
density- independent regime at an observer time t = 10 mins. In addition, we know that 
(l+Zc) OC t^^u^"^ based on eqn (12). As the observer time or observer's frequency is increased, 
the critical redshift decreases, if the other parameters remain unchanged. Therefore, in the 
range of redshifts concerned, GRBs in these two different density profiles are always in the 
density-independent regime and the corresponding fiuxes will always be same. On Figure 
5 (a) and (b), the X-ray flux is calculated for observer times t — 8 hours, 12 hours, 1 day 
and 2 days for Chandra and XMM, respectively. The flux is integrated over the 0.4 — 6 keV 
range of the Chandra ACIS instrument, and over 0.15-15 kev for XMM. Although it takes 
about 1 day or so for Chandra to slew onto the source, we see that Chandra is still able 
to detect GRBs with the typical parameters considered here up to z ^ 30 at 1 day with a 
10 ks integration. XMM has proved itself capable of < 8 hours slewing onto GRBs and 
has a similar sensitivity to Chandra, hence it might be able to detect higher redshift GRBs 
than Chandra does. In Figure 5(c), we have integrated the X-ray flux over the Swift XRT 
frequency range of 0.2 — 10 keV. Although Swift XRF has a relatively lower sensitivity than 
Chandra and XMM, this is compensated by its quick slewing time, less than 1 minute. In 
Figure 5(c) wc sec that Swift XRT can easily detect typical GRBs up to redshifts ~ 30, 
if they are observed within 1 hour after the trigger. Therefore, if GRBs exist at very high 
redshift, we can expect these detectors to be able to measure them in x-ray band. 



4. Summary and Conclusions 

In this paper, we have calculated the spectral time evolution and the flux in the near- 

IR K- and M-bands as well as in the X-ray band from GRBs at very high redshifts and 
different times. In previous work, Ciardi & Loeb (2000) calculated IR fluxes as a function 
of redshift and time for the standard forward shock model of afterglows. Lamb & Reichart 
(2000) discussed the optical/lR as well as gamma-ray fluxes of some observed bursts at an 
observer time of 1 day when placed at different redshifts. Here we have introduced several 
new elements in our analysis, motivated by recent developments in the observations as well 
as in the modeling of bursts. The most important of these is that we consider, in addition to 
the forward shock, also the reverse shock, which has now been inferred in three GRBs from 
prompt follow-ups. The quick response capability of a number of ground and space observing 
facilities coming on-line in the near future means that one is far likelier to observe the early 
stages of the GRB and of its afterglow evolution. Thus there are excellent prospects for 
observing the reverse shock thought to be responsible for the prompt optical flashes, which 
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is prominent only at early times in the burst evolution. The observation of reverse shocks 
will, in addition, provide significant independent information on early GRB evolution, such 
as the initial Lorenz factor, the strength of magnetic fields, etc. (e.g. Zhang et al. 2003). 
Another difference with previous fiux calculations is that we have used significantly updated 
model parameters, based on new data acquired in the past two years. Thus, for instance, we 
make use of the emerging consensus view that esj is usually smaller than eg,/. As is seen 
from the calculations presented here, the magnetic field equipartition value has a significant 
effect on the flux. 

We have calculated the fiux from high redshift GRBs taking typical parameters, which 
gives a sense for how far GRBs can be detected using current or forthcoming instruments. 
In reality, these parameters would vary over a wide range, which would also affect the 
detectability. Zhang et al. (2003) have reexamined the three cases of GRB 990123, GRB 
021004 and GRB 021211 and found evidence for an enhancement of the magnetic field in the 
reverse shock over that in the forward shock, by as much as a factor of TZb = Br/Bf f« 15 
(GRB 990123), i.e. sb in reverse shock is much larger that in the forward shock. If so, based 
on eqns. (A12-A17), we can see that if TZb increases, the observed (optical/IR) fiux for the 
reverse shock component increases significantly as illustrated by comparing Figures (2) and 
(3). Separately, it has recently been recognized that there are burst to burst variations in 
the total beaming-corrected energy, and thus, the forward shock peak afterglow fluxes may 
also be significant (Bloom, Frail, & Kulkarni 2003). The so-called 'fGRBs' exhibit rapid, 
jet-break like decays at early times, before the < 1 day point at which they would usually be 
expected. Thus we can expect that the reverse shocks in this kind of GRBs are less bright. 

Most of the radiation in the optical and ultraviolet bands from high redshift extragalactic 
sources, including GRBs, is absorbed by the intergalactic medium and the diffuse gas in our 
own galaxy. The X-ray and infrared bands are therefore of major importance for detecting 
and tracking high-redshift GRBs. Several major ground-based telescopes as well as smaller 
robotic facilities have or will have infrared sensitivity in the K, L and M-bands. Next 
generation spacecraft such as Swift have X-ray and optical/UV detectors, while the James 
Webb Space Telescope (JWST) frequency range extends out to 27 /xm, being most sensitive 
in the 1-5 /xm J,H,K,L and M bands. In the X-ray band, the Chandra and XMM sensitivities 
in 0.2-10 keV are substantially higher than that of Swift XRT, but their slewing time (< 1 
day) limitations make Swift XRT a unique instrument for X-ray follow-up during the first 
day after a GRB trigger, when the burst is brighter. In spite of this, all three spacecraft 
should be able to detect very distant GRBs, if they exist, e.g. at z > 30. 

For the nominal GRBs considered here, the luminosities are comparable to those of the 
currently detected ones. According to theoretical modeling, the fractional number of GRBs 
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expected at 2; > 5 is > 50% of which ~ 15% may be detectable in flux-hmited surveys such as 
Swift's (Bromm & Loeb 2002). It was reported that HETE-2 sees 13 out of 14 GRB optical 
afterglows by Sep. 2003, which means that the high-z GRB fraction is small (Ricker 2003). 
However, recently another high-z GRB candidate (GRB 031026) was detected and proposed 
(e.g. Attcia et al. 2003), which increases the high- 2; fraction to be close to 15%. We note 
also that at 2; > 6 — 10 the first generation of (pop III) stars are likely to lead to black holes 
with masses 10 — 30 Mq, and hence to GRBs whose luminosities could be factors 10-30 times 
higher (e.g. Meszaros & Rees 2003) than assumed here. In this case, the fraction at 2; > 5 
detectable in Swift's flux-hmited survey could be > 20 — 30% of the total, or > 20/year. 

In the K and M bands (2.2 and 4.8 /xm) the JWST and other telescopes should be able 
to detect afterglows out to 2; < 16 and 33 within observer times 1 day for integrations times 
(with JWST) of 1 hour (at a resolution R=1000 and S/N=10; see Figure 3 and Figure 4). 
These bands are accessible also to ground-based telescopes already on line, before the JWST 
launch. The effect of reverse shocks, which are brighter in the 0/IR at early source times 
for some bursts, makes for a significantly increased sensitivity at high redshifts at observer 
times < 1 day. 

We have also considered the effect of two different types of the near- burst environment, 
one assuming that the external density evolves with redshift similarly to that in protogalactic 
disks, and the other assuming approximately redshift-indcpcndcnt conditions regulated by 
radiation pressure, based on primordial star formation calculations. The predicted X-ray 
fluxes, being due mainly to forward shocks above the cooling frequency, arc independent of 
the external density regime, as well as insensitive to the existence of reverse shocks. However 
the IR fluxes are sensitive to which of these density regimes prevails, and at early times are 
very sensitive to the presence and strength of a reverse shock component, in particular at 
early times and redshifts ^15. Combining these two types of IR and X-ray flux information 
will thus provide very important tools for detecting GRBs (if present) out to very high 
redshifts, for studying their local environments, and for investigating the effects of reverse 
shocks as well as the prompt phases of the bursts and their afterglows. 

We are grateful to the referee for useful comments, to Shiho Kobayashi, George Chartas, 
David Burrows, Jian Ge and Zheng Zheng for useful discussions, and to NASA NAG5-9153, 
NAG5-9192, NAG5-13286 grants and the MoneU Foundation for support. 
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A. Appendix: Flux Evolution 

A.l. Forward Shock 

The radiation emitted by a source at a redshift z at frequency Us over a time 5t will 
be observed at z=0 at a frequency ^(^{z = 0) = ^s/{^ + over a time 5tQ{z = 0) = 
(1 + z)5ts. The luminosity distance for a flat universe + Q,m — 1, = 0.27 and 
Hubble constant (i7/100) km/ s/ Aipc = O.lhjo can be approximated, in units of 10^^ cm, as 
D28iz) ^ 4.49(1 + z)[l- 1.115(1 + z)-'^/'^] (Pen 1999). Substituting this redshift dependence 
and eqn. (5) into equations (1) and (2), we have for the fast cooling case 

r CesjEli't'/V^^iyy^l + z)-V6^-2 ^ < ^^^^ 

= \ C,e-fEli%'/'u-y\l + z)-'/^x~' u^j < v < u^j (Al) 

[ C3eg'//^er^)£;(^+')/'i-(=^^-')/V-^/2(l + z)^-^)l'x-' v^j < v 

For the slow cooling case, we have 

= < C5ei"/'^/'e(^-^)E(r^)/^t-='(^-^)/^nVV(i-f)/2(l + z)^p-^)''x-^ i^mf<i^<i^cf 
[ C3eg'//^ei^-^)£;5^+')/^t-(^^-')/^V-^'/^(l + z)(P-'y^x-' u,f<u 

(A2) 

where Ci=0.402, C2 = 8.65 x 10^ C3 = (8.6 x 10^) x (4.6 x io14)(p-i)/2 = § 5 x lO^o (p = 2.5), 
C4 = 0.071,and C5 = 5.5 x 10^ x (4.58 x iQ^^)(p-^)/2 = 5.4 x 10^^ (p = 2.5) and x is defined 
as 

a; = [1-1.115(1 + z)-^/2]. (A3) 

Here, esj, Ce, -E'52 have been assumed to be constant parameters, while others like D2s{z) and 
n may be redshift dependent in different cases, u is the observed frequency. For simplicity, 
we present here the scaling relation for the flux with several parameters which may change 
with redshift. We substitute each of these quantities into the expressions above, and get the 
scaling relations with redshift for different cases. 

Fast cooling case: 

r n5/6(l + z)-5/6a;-2 

oc I n0(l + ^)-5/4a;-2 

[ nO(l + ;2)(P-6)/^x-2 



UcJ <U < UmJ (A4) 



Slow cooling case: 
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r nl/2(l + ^)-V6a;-2 

oc I ni/2(i + ^)(f-5)/4^-2 
i n0(l + ;2)(^'-6)/4x-2 

Substituting the redshift dependence of the number density into the relations above 
gives straightforwardly the scahng relations for the different density cases. 

A. 2. Reverse Shock 

The scalings here are taken from Kobayashi (2000) and Zhang et al. (2003). Here 
the crossing time is defined as = '>Tiax{T,tdec), where T is the burst duration and td^c is 
defined as [{3E /Anrj'^nmpC^y/^ /2r]'^c]{l + z) in the observer frame. Also, rj — mm{r],r]'^/r]), 
where r] is the initial Lorentz factor and r]c is defined as the critical initial Lorentz factor 
77c = UbEli^n-^/^T'^^^i^f/^ (Zhang et al 2003). For the thin shell case, one has = W 
and "fj — T), while for the thick shell case, one has tx —T and rj — t}c/v- 

In the thick shell case, the typical parameters at crossing time tx are: 

where TZs = Br/Bj = {eB^r/ (^Bjf^ ■ 

The scaling relations before and after the shock crossing time tx are 

t <i tx '■ ^m,r t , I'cr t , Fj^ffy^ j. (X t ^ (-'^''^) 



t>tx: Vm,r OC r^3/48 ^ ^-3/2^ ^^^^ ^ ^-73/48 ^ ^-3/2^ p^^^^^ ^ ^-47/48 ^ ^-1 (^8) 

For the thin shell case, at the crossing time, one has, 

i^m,r{tx) ^V'^^BJ^mjitx), i^c,r-(^x) = ^B~V,/(^x), F^^rnAtx) ^ r]T^BFu,m,f{tx) (A9) 



I'mj <y < i^cf (A5) 



The scaling relations before and after the shock crossing time tx are 
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t<ty,: Vm,r = ^^ ^c,r « t ^ F^,^,^ OC t^/'^ (AlO) 
t > : ^n.,. OC t-''/'' ^ t-^l\ V^, OC r^V35 ^ ^-3/2^ ^ ^-34/35 ^ ^-1 (^11) 

Before the crossing time, for the thin shell case observed flux can be expressed as: 



\F,,iv]J^v-'l\fiRs-^'\tlUfl''\ V,, <v< (A12) 



'va,r 



while for the thick shell the flux is: 



c,r 



V <v. 



m.r 



= { [F.,f^i^7'^^'^-^-'^/'][^^'-^^^B(^+^^/'(iAx)^/^] ym,r < y < yc,r (Al5) 

[Fu.y^j'^^'y'Sy-'/'] [e-'^Us^"-'^/'] u,,r < y 

After the crossing time, for both the thick shell and thin shell cases the expressions for 
the flux are: 

[F.,ty;J^'y'/']mB'{t/t,)-'/'] y < yc,r 
F^= { yc,r<y< ym,r (A16) 

ym,r < y 
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m,r 



[^.,f^i7'^/'^-^-'^/'][^^'-^^^B(^+^^/^(V^x)(^-=^^)/1 ym,r < y < yc,r (Al7) 

i/^ . < V 



Substituting the redshift dependence into equations above for the reverse shock, we 
obtain similar scahng relations as those for the forward shock. 

For the thin shell case before crossing time: 

Fast coohng case: 

cx <j ni3/i2(i + ^)-5/4a.-2 ^^^^ < < ^^^^ (^18) 



Slow cooling case: 



n 



m.r 



OC { n(3P+3)/4(l + ^)(p-5)/4^-2 ^^^^ ^ ^ ^ ^^^^ ^^19) 

^(5p-4)/4(^^^)(p-6)/4^-2 ^^^^^^ 

For the thin shell case after the crossing time: 
Fast cooling case: 



nlVl8(l + ^)-5/6^-2 

F^ (X { U^r <y < ym,r (A20) 

, Ur, 

Slow cooling case: 

r nV6(l + ^)-7/6^-2 l/<i/^,, 

OC I nV3(i + i/^,^ < u< Uc,r (A21) 

I Uc,r < y 

One sees that if the number density around the GRBs does not change with the redshift, 
i.e. n — constant for all z, the reverse shock emission depends on redshift in the same way as 
the forward shock emission. However, in the n oc (1 + 2;)^ case, the behavior for the reverse 
and forward shock emission is different. The thick shell behavior can be obtained in a similar 
way. 
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Fig. 1. — Typical light curves, for a rcdshift z = 1. Reverse shock emission (dashed), forward 
shock emission (solid), total flux (symbols). Parametersie^j = 0.001, 7^^ = Br/Bf = 5, eg = 
0.1,^52 = lO,p = 2.5,7]^ 120, no = 1 cm'^. a): V band {u = 5.45 x 10^^ Hz); b): K band 
(i/ = 1.36 X 10^^ Hz). 
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Fig. 2. — Combined forward and reverse shock observed flux as a function of redshift for 
^Bj = 0.025 and TZb = B^/Bj = 1. Forward shock (symbols), reverse shock (without 
symbols). Solid, dashed and dotted lines indicate emission at different observer times t=10 
mins, t=2 hour and t=l day respectively, a): V-band {u — 5.45 x 10^^ Hz); b): K-band 
{u = 1.36 X 10^^ Hz); c): M-band {u = 6.3 x 10^^ Hz). Straight lines: In V band sensitivities 
are for ROTSE at very early and late times; in K and M bands sensitivities for JWST K 
& M bands are estimated for a resolution R=1000, S/N=10 and integration time of 1 hour. 
Parameters: n — 1 cm~^, — 0.1, £^52 = 10, p = 2.5, r] = 120. 
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Fig. 3. — Combined forward and reverse shock observed flux as a function of redshift for the 
two density profiles n = Uq = constant ( without symbols) and n = no(l + z)^ (symbols) 
with Uq — 1 cm~^. solid, dashed and dotted lines indicate emission at different observer 
times t=10 mins, t—2 hour and t=l day respectively, a): V-band {u — 5.45 x 10^^ Hz); b): 
K-band {u = 1.36 x 10^^ Hz); c): M-band {u = 6.3 x 10^^ Hz). The limiting ROTSE and 
JWST sensitivities are the same as in Figure 2. Parameters: cbj — 0.001, 7?.^ = B^/Bf — 
5,ee = 0.1,£;52 = 10,p = 2.5,77 = 120. 
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Fig. 4. — Combined forward and reverse shock observed flux as a function of redshift for 
the two density profiles n = tiq = constant (without symbols) and n = no(l + z)'^ (symbols) 
with no — 0.01 cm"^. Sohd, dashed and dotted lines indicate emission at different observer 
times t=10 mins, t—2 hour and t=l day respectively, a): V-band {u — 5.45 x 10^^ Hz); b): 
K-band {u = 1.36 x 10^^ Hz); c): M-band {u = 6.3 x 10^^ Hz). The limiting ROTSE and 
JWST sensitivities are the same as in Figure 2. Parameters: cbj — 0.001, 7?.^ = B^/Bf — 
5,ee = 0.1,£;52 = 10,p= 2.5,7y = 120. 
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Chandra in 3.6 ks & 10 ks integration time@ 0.15-15 Kev 



XMM in 5 ks & 10 ks integration time@ 0.15-15 Kev 
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Swift in 5 mins & 1 ks integration time@ 0.2-10 Kev 
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Fig. 5. — Observed X-ray fluxes for GRB afterglows at different redsiiifts, integrated over 
the observing energy ranges of 0.4-6 keV for Chandra, 0.15-15 Kev for XMM and 0.2-10 
keV for Swift, respectively. The emission is in the density-independent regime, above v^j. 
a): for Chandra, the fluxes for the observer times tohs — 8 hour, 12 hours, 1 day and 2 
days as compared to its sensitivities shown as horizontal lines for integration times of 3.6 ks 
(dashed) and 10 ks (solid), b): for XMM, same observer times as Chandra's. The sensitivity 
horizontal lines are for integration times of 5 ks (dashed) and 10 ks (solid), c): for Swift 
XRT, the fluxes are for observer time t^^^ = 10 mins, 20 mins, 1 hour. The sensitivity 
horizontal lines are for integration times 300 s (dashed) and 1 ks (solid). Parameters: esj = 
0.001, = Br/Bf = 5,ee = 0.1,^52 = 10,p = 2.5,77 = 120. 



